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We study gravitational lensing by clusters of galaxies in the context of the generic 
class of unconventional gravity theories which describe gravity in terms of a metric and 
one or more scalar fields (called here scalar-tensor theories). We conclude that if the 
scalar fields have positive energy, then whatever their dynamics, the bending of light 
by a weakly gravitating system, like a galaxy or a cluster of galaxies, cannot exceed the 
bending predicted by general relativity for the mass of visible and hitherto undetected 
matter (but excluding the scalar field's energy). Thus use of general relativity to in- 
terpret gravitational lensing observations can only underestimate the mass present in 
stars, gas and dark matter. The same conclusion obtains within general relativity if a 
nonnegligible part of the mass in clusters is in the form of coherent scalar fields, i.e. 
Higgs fields. The popular observational claim that clusters of galaxies deflect light much 
more strongly than would be expected from the observable matter contained by them, if 
it survives, cannot be interpreted in terms of some scalar-tensor unconventional gravity 
theory with no dark matter. And if the observations eventually show that the matter 
distribution inferred via general relativity from the lensing is very much like that de- 
termined from the dynamics of test objects, then scalar-tensor unconventional gravity 
will be irrelevant for understanding the mass discrepancy in clusters. However, even a 
single system in which the dynamical mass determined from virial methods signiflcantly 
exceeds the lensing mass as determined by general relativity, would be very problem- 
atic for the dark matter picture, but would be entirely consistent with unconventional 
scalar-tensor gravity theory. 
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1. INTRODUCTION 



The matter content of extragalactic systems, e.g., galaxies and clusters of galaxies, 
can be assessed in two independent ways: from the dynamics of test objects, e.g., 
extended rotation curves of disk galaxies delineated by neutral hydrogen, or the velocity 
dispersion profile of cluster delineated by galaxies, and from the deflection and focusing 
of electromagnetic radiation, e.g., gravitationally lensing cluster of galaxies. Because the 
"bending of light" by gravitational fields is intrinsically a relativistic effect, the second 
approach provides a way to test the relativistic aspects of gravitation at the extragalactic 
level. This is particularly important in light of a variety of proposals (Tohline 1982, 
Milgrom 1983a,b,c, Bekenstein & Milgrom 1984, Sanders 1984, Bekenstein 1987, Kuhn 
& Kruglyak 1987, Sanders 1986a,b, 1988, 1989, Manneheim & Kazanas 1989, Kazanas 
& Manneheim 1991, Sanders 1991, Bekenstein 1992, Romatka 1992, Milgrom 1993, 
Kenmoku et al. 1993) that seek to replace the dark matter paradigm with one or 
another deviation of gravity from standard theory at extragalactic scales. Whatever 
the success of such schemes, and there have been many in the case of Milgrom's, they 
must still face the challenge of providing a theory of light bending that correctly predicts 
the features of gravitational lenses, arcs, etc. 

Now in the standard gravitational theory (general relativity [GR] and its Newtonian 
limit), the two approaches to measure the matter content of a system use the same tool 
to describe the predictions: the Newtonian gravitational field g. Thus the acceleration 
of a test object in a galaxy or in a cluster of galaxies is just this g evaluated at the 
position of the object. And the bending angle of a light ray as it passes near the system 
along an approximately straight path may be stated to lowest order in where v 

denotes the typical velocity in the system, as 



where ± denotes the component normal to the ray's direction, and dz is the element 
of length along the ray. In this context do observations of cosmic gravitational lenses 
affirm the dark matter paradigm in the context of standard gravity theory and rule out 
novel gravitational theories as alternatives to dark matter ? 

Unfortunately, up to now both observational and theoretical ambiguities have stood 
in the way of a sharp answer to this question. As detailed below, the observations of 
lensing systems, principally clusters of galaxies, do not yet permit a very accurate com- 
parison of the lensing mass with the standard dynamical mass derived from kinematical 
measurements. And unconventional gravity theories have not, up to now, led to model 
independent predictions for the lensing that can be confronted with the evidence and 
the predictions from the dark matter picture. For instance, it has not been known in 
any generality whether in scalar-tensor (ST) type gravitational theories without dark 
matter, the scalar field is expected to bend light beyond what is accomplished by the 
luminous matter. We now review the status of the observations and unconventional 
gravity theoretical frameworks. 

Regarding observations of gravitational lenses, in those cases of strong imaging 
events i.e., the formation of multiple images or Einstein rings, where the lens is an 
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individual galaxy (see Blandford & Naryan 1992), there is no indication of a discrep- 
ancy between the lensing and luminous masses (Breimer & Sanders 1993). This is not 
surprising because, typically, the impact parameter of photons forming multiple images 
is only a few Kpc, and there is convincing evidence that the bright inner regions of 
galaxies are typically dominated by the visible component (Kent 1987, van Albada & 
Sancisi 1987). True, Maoz and Rix (1993) have claimed that the observed frequency 
and image separation of multiply-lensed quasars can only be achieved if galaxies - in 
particular elliptical galaxies - have dark halos i.e., the dark component makes a contri- 
bution to the image magnification and separation. However, this result depends upon 
rather uncertain model parameters and scaling relations for elliptical galaxies. 

The issue of lensing by dark haloes has also been examined by observations of the 
lensing effects of forground galaxies on background galaxies within projected separations 
between 10 kpc and 200 kpc i.e., weak gravitational imaging (Tyson et al. 1984). In 
this situation of large projected separations, the presence of the foreground mass does 
not create multiple images of a background galaxies, but does systematically distort 
their shapes. This is very relevant to the problem of the mass discrepancy because it is 
a light-bending probe of the mass distribution in the outer parts of galaxies where the 
dark halo is thought to dominate. In the present context it is interesting that the results 
of Tyson et al. are consistent with less extensive dark haloes than implied by measures of 
the extended rotation curves. For example, these results imply that a typical lens galaxy 
has a total mass less than 3 x lO"*^-*^ Mq within a radius of 80 kpc; this is comparable to the 
luminous mass of a typical elliptical galaxy i.e. a galaxy with the fiducial luminosity, 
L*, in the Schechter (1976) luminosity function. Is this result a falsification of dark 
matter hypothesis ? Unfortunately the Tyson et al. result remains controversial. For 
example, Kovner & Milgrom (1987) argue that the assumption of infinite distance for 
the background sources underestimates the total mass of the forground galaxies by a 
factor of three. On the other hand, Breimer (1993), using a realistic distribution of 
redshifts for the background sources, has confirmed the result of Tyson et al. provided 
that the intrinsic source size is larger than the seeing disk. Clearly, further observations 
of this sort under conditions of good seeing would be most useful. But it can be said 
that, up to now, there is no conclusive evidence for lensing by dark matter in individual 
galaxies. 

The most persuasive evidence in favor of lensing by dark matter is provided by 
the luminous arcs seen in the central regions of a dozen rich clusters (Blandford & 
Naryan 1992). It is now well-established that these are background galaxies imaged 
into partial Einstein rings by the forground cluster (Lynds & Petrosian 1986, Soucail 
et al. 1987). The presence of one of these arcs implies, via eq. (1), the existence 
of more than 10^"^ Mq projected within the inner 200 kpc (the typical Einstein ring 
radius), masses which exceed, by a factor of ten or twenty, the projected luminous mass 
(Bergmann, Petrosian & Lynds 1990). Moreover, the observations of systematically 
distorted images of background galaxies within projected radii of 1/2 to one Mpc in 
two distant clusters (Tyson, Valdez, & Wenk, 1990) apparently requires a substantial 
extended dark component. In the cluster A 370, the most well-studied of the arc 



systems, the cluster velocity dispersion predicted by the observed arc is consistent within 
the observational errors with the observed velocity dispersion of the cluster galaxies, 
assuming that the visible galaxies trace the dark mass distribution and that the velocity 
distribution is isotropic (Mellier et al. 1988, Breimer & Sanders 1992). There thus 
appears to be consistency between the virial mass and the lensing mass. This would 
seem to be evidence in favor of the dark matter hypothesis. 

However, these observations do not yet rule out generic unconventional gravity the- 
ories without dark matter for two reasons. First of all, determinations of the mass 
and mass distribution by standard dynamical methods are notoriously crude. Virial 
methods, including those based upon the Jeans equation, which use measurements of 
the line-of-sight velocities of cluster galaxies are beset by uncertainties resulting from 
limited statistics, the unknown degree of anisotropy in the velocity distribution, the the 
effects of subclustering and contamination, and the unknown distribution of the total 
(observed and unobserved) mass (The and White 1986, Merritt 1987). The potentially 
powerful method of using the observed density and temperature distributions of the hot 
X-ray emitting gas has been limited by the lack of high resolution spectral information; 
i.e., the temperature distribution is very uncertain (Sarazin 1988). For these reasons, 
dynamical masses are typically determined to an accuracy of a factor of two at best. 
Secondly, recent reanalyses of the gaseous component of clusters on the basis of results 
from ROSAT and earlier X-ray observatories indicate that the contribution of the hot 
gas to the total cluster mass may be greater than previously supposed - up to 40% in 
some cases (Hughes 1989, Eayles et al. 1991, Briel, Henry & Bohringer 1992, Bohringer, 
Schwarz & Briel 1993). This goes a long way towards relieving the discrepancy in clus- 
ters. Moreover, many rich clusters show evidence for cooling flows which can deposit up 
to 10''"'^ Mq in the central regions in a Hubble time, presumably in the form of cold gas 
or low mass stars (Sarazin 1988). This raises the possibility that, in clusters containing 
giant arcs, the lensing mass is nothing more exotic than hot gas or objects formed out 
of the hot gas. Future X-ray satellites, which will probe the temperature distribution 
of the gas in clusters, should clarify the issue of the mass distribution in the central 
regions and the contribution of the hot gas to the total mass in the lensing clusters. 

In light of the above discussion, claims sometimes heard, e.g. Dar (1993), that the 
extant observations of gravitational lenses, when compared with dynamical determi- 
nations of the virial mass, certify the validity of general relativity in the extragalactic 
regime are evidently premature. They do not reckon with the intrinsically crude determi- 
nation of g provided by the dynamical measurements. Neither do they take cognizance 
of the dearth, in the extant literature, of clearcut predictions for gravitational light 
bending in theories of gravity, which aim at explaining the data without an appeal to 
dark matter. There is a good reason for this dearth. Bending of light in a gravitational 
theory depends on details of the theory: what fields constitute gravity, how are they 
coupled to gravitating matter, what equations are satisfied by these fields, . . . ? On the 
face of it, predictions of light bending would have to be worked out separately for each 
theory, and it certainly could not be ruled out a priori that some model theory could 
produce light bending of the same strength as GR. 



The actual situation is not so bad. In many known gravity theories where a scalar 
field plays a gravitational role, it does so by entering as a conformal factor in the relation 
between the primitive (Einstein) metric and the physical one. This fact implies, in a 
rough way, that the scalar field will not induce any light bending of itself (Bekenstein 
1992, Romatka 1992). Two advances were required in order to turn this observation 
into the solid conclusion (this paper) that ST unconventional gravity theories cannot 
reconcile the lensing and dynamical (velocity dispersions) data for clusters without 
requiring much ordinary or exotic unobserved matter. 

One advance was to establish the generic way to couple a scalar field to the metric. 
Bekenstein (1993) has shown that the relation between Einstein and physical metric 
may be generalized to the form of a disformal transformation in which the stretch of 
length in the special spacetime direction delineated by the scalar field gradient is either 
larger or smaller than the average one. Because of the anisotropic stretch, it is no longer 
automatic that the scalar field cannot cause bending of light. The disformal transfor- 
mation, discussed in Sec. 2., seems to be the generic way to introduce a scalar field into 
a gravity theory which has a fighting chance to pass the tests of GR. Considerations 
of causality dictate (Bekenstein 1992) that the special direction undergoes a smaller 
stretch than the average over all directions. This particular sign turns out to be critical 
in the results to be described presently. 

The second advance required was a relativistic calculation of light bending for a 
generic ST gravitational theory. The rough insight that, as a conformal factor in the 
physical metric, a scalar field cannot bend light, does not reckon with the fact that the 
stress energy of the scalar field can, by way of its contribution to the sources of the 
Einstein metric, contribute to bending of light. This contribution cannot be ignored: 
the evidence is that there is a significant non-luminous mass in clusters, which could 
very well turn out to be energy of coherent fields. A truly relativistic calculation of light 
bending will take into account this cff'ect. Such is carried out in this paper, in Sees. 3. 
and 4. for the bending due to a spherically symmetric system in a generic gravitational 
theory involving scalar fields with positive energy densities. 

The result (Sec. 5.) is surprising: despite all the subtleties mentioned, the bending 
angle cannot exceed that predicted by the standard theory for the actual matter (ex- 
cluding the scalar field) present in the system. Thus it is a theorem that in a generic 
ST theory of gravity, the scalar field cannot enhance lensing. However, if a theory is 
presented as a candidate for gravity with no dark matter, it will, almost by construc- 
tion, predict that the visible matter in the system generates a stronger gravitational 
field than would be predicted by standard theory. Therefore, modified gravitational 
theories based on metric and scalar fields and assuming no dark matter, and standard 
theory with dark matter must give distinctly different predictions for the comparison 
of dynamical and light bending data. That is, a modified gravity theory of the type we 
have in mind cannot mimic all predictions of standard theory with dark matter. 

Thus, as we conclude in Sec. 6., if future measurements of lensing and of velocity 
dispersions in clusters (or of the temperature distribution of the hot gas) certify the 
presently popular view that gravitational lensing requires the same mass discrepancies 



in galaxy clusters as are implied by test body dynamical measurements, then unobserved 
matter makes up a large fraction of the mass in such clusters, whether gravity is best 
described by conventional theory or by some unconventional ST theory. In the later 
case the indication would be that in the cluster regime the predictions of the unconven- 
tional theory differ little from those of conventional gravity theory, with the difference 
between the theories being significant, if at all, only on scale of galaxies. By contrast, 
if future observations show that the lensing by clusters suggest masses below those im- 
plied by the velocity dispersions of galaxies or the distribution of hot gas, then dark 
matter would loose much of its appeal as an explanation for the entire mass discrepancy, 
while unconventional gravity theories of ST type would become quite relevant. In this 
light gravitational lensing is seen to be a crucial tool for the issue of dark matter vs. 
unconventional gravity. 

In what follows we choose units such that the speed of light c = 1, but we shall retain 
Newton's constant explicitly. Our conventions for gravitational theory are those of Mis- 
ner, Thorne and Wheeler (1973). In particular we assume metric signature (— , +, +, +). 

2. CHARACTERIZING GRAVITATIONAL THEORIES 

2.1. From General Relativity to Scalar-Tensor Theory 

Since we shall discuss gravitational theories which involve a scalar field, it is in order 
to start with a discussion of why a scalar field is the most natural entity that can be 
used to modify standard gravitational theory. The full characterization of a gravitational 
theory must be relativistic. The simplest such theory, GR, is formulated in terms of 
the metric of spacetime, g^^p, which enters into all dynamical equations for matter 
in gravitational fields in the same manner that Minkowski's metric enters in gravity's 
absence. The simplicity of GR is manifested in the absence of other gravitational fields 
apart from the metric, and in the fact that the dynamical action for the metric is the 
simplest invariant action that can be built in 4-D curved spacetime: 

= [ Ryf^dt^X (2) 

Here R is the scalar curvature built from the metric ^Tq,^ and its first and second deriva- 
tives, and g is the metric's determinant. 

The success of GR in confronting the results of solar system precision experiments 
(Will 1992) suggests that the correct gravitational theory, if not identical to GR, must 
nevertheless share with it some key features. Therefore, any theory which attempts to 
displace dark matter as an explanation of the mass discrepancy in large astrophysical 
systems must show clear kinship with GR. The action (2) and the underlying metric 
Qo^jj are such key features of GR that it is tempting to see them as required of any 
competing theory. We shall assume that a theory of gravity, to be credible, must have 
action (2) in some choice of local units (in some conformal frame). 

Similarly, the role of a curved spacetime metric in rewriting special relativistic 
physics to include the effects of gravity seems a required feature of theories that are 
to compete succesfuUy with GR. We shall assume, therefore, that the generic theory to 



be considered uses some such metric. However, it does not follow that this last metric, 
call it gajji is identical to g^y^p. In GR that is the case, and the equivalence is a facet 
of the strong equivalence principle which is part and parcel of GR. However, it is only 
the Einstein equivalence principle, a weaker principle, which is supported by many ex- 
periments. As stressed by Dicke long ago (Dicke 1962), those experiments are evidence 
only for the existence of the metric ^q,^, and it is logically consistent with all extant 
experiments for the metrics ^q,^ and ^Tq,^ to be distinct. The metric g^^p^ often referred 
to as the Einstein metric, shall here be called the gravitational metric; ^q,^ shall here 
be referred to as the physical metric. This last is the metric determined by rods and 
clocks built of matter. 

In the ancient Nordstrom theory of gravity (Nordstrom 1913), which actually pre- 
ceded GR, the metric ^Fq,^ was taken as flat, whereas ^q,^ was taken as conformal to it, 
with the conformal factor being the square of a scalar field which obeyed the standard 
massless scalar equation written with ^Tq,^. Thus Nordstrom's theory is a pure scalar 
gravity theory. Despite its inherent simplicity, the theory was rejected early. It is amus- 
ing, given our concerns here, that its failing was predicting that there is no bending of 
light. 

The Bergmann-Wagoner scalar tensor gravitational theories (Will 1992), of which 
Brans-Dicke theory is the simplest, are a natural marriage of Nordstrom's theory with 
GR. The gravitational metric g^^j^ is taken to satisfy Einstein-like equations, whereas the 
physical metric g^^^ is taken to be conformally related to it, with the conformal factor 
being intimately related to the dynamical scalar field of the theory, which together with 
the metric constitutes gravity. 

2.2. Disformal Transformation 

Of course, once a scalar field, if), is admitted into gravitational theory, the relation 
between physical and gravitational metrics may be more general. For example, consider 
a disformal relation (Bekenstein 1992) between the metrics, 

-g^p = exp(2V') [A{I)g^p + L^B{I)i^,o^ V',/3 ] (3) 
where L is a scale of length, and A and B are two functions of the invariant 

/ = LV^V,aV',/3 (4) 

The form of the factor exp(2?/') is dictated by the reasonable requirement that a shift 
in the zero of •0 shall have no physical consequences. Indeed, adding a constant to '0 
merely changes the global units of length paced out by the physical metric. 

The relation (3) is the most general between g^j^ and gg^p based on a scalar field 
if) which involves only the first derivatives of ij). Generalizations of the relation (3) 
involving second derivatives of '4> obviously entail a higher derivative theory, i. e. one in 
which the matter field equations involve field derivatives of third or higher order. It is 
well known that such theories tend to display causality problems such as preacceleration 
and runaway solutions, and it seems wise to exclude them at the outset. It has also 
been shown (Bekenstein 1993) that the disformal relation (3) is the most general relation 



between a given gravitational metric ^q,^ and a Finsler geometry for the matter equations 
which is aUowed by the requirements of causahty and universahty of free faU. 

Let us recapitulate the structure of the generic ST gravitational theory. The gravi- 
tational action (2) is responsible for dynamics of the gravitational metric gfQ,^. The full 
Einstein-like equations are obtained from the variation of the full action 

'S'tot = Sg + Sm + (5) 

with respect to g^.^. Here S^n is the action of the matter fields obtained by replacing 
the Minkowski metric in the special relativistic action for matter by ^q,^, and partial 
derivatives by covariant derivatives constructed with ^q,^. And is an action for if; 
built out of -0, its first derivatives, and g^p. The dynamics of ifj comes from its field 
equation which is obtained by varying S'tot with respect to •0. They combine with the 
dynamics of g^.^ coming from the Einstein-like equations, to determine how g^,^ varies 
in spacetime. 

2.3. Examples and Generic Theory 

Let us look at some examples of theories encompassed by the stated framework. 
First, GR has ~ with '0 = by fiat, with A = B = 1. Here the two metrics are 
equivalent. Brans-Dicke theory has A = 1, B = and 

= J {00 + 3/2)g^^i;,a i^,pV^d^x (6) 

with A = being the scalar field in Dicke's (nonmetric) form of the theory (Dicke 1962), 
and u> being the usual Brans-Dicke parameter (really an inverse coupling constant in 
modern jargon). As is well known, Brans-Dicke theory is in good accord with all solar 
system tests of gravitational theory for oj > 10^. The Bergmann- Wagoner theories differ 
from the special Brans-Dicke case in that oj = (jj{ip) and a potential term is added to 
the kinetic part of the action (6). In the extragalactic context all the above theories are 
esentially indistinguishable. 

Several ST theories have been proposed as alternatives to the dark matter hypothesis 
coupled to standard gravity theory. First in simplicity is Sanders' two scalar field theory 
(Sanders 1986), which implicitly assumes A = 1 and 5 = 0, and in which a linear 
massless scalar field is accompanied by a massive scalar field whose kinetic action has a 
sign opposite to the usual one. This negative energy field is repulsive. The theory may 
pass the various tests of GR, and is not obviously in contradiction with the universality 
of free fall. Yet, the negative energy is a detractive feature. 

Following Sander's theory in complexity is the aquadratic lagrangian (AQUAL) 
theory (Bekenstein and Milgrom 1984). It has A = 1, B = and 

= / F{I)^gd^x (7) 

with F{I) oc fi/"^ for small / and F (x I for large /. This theory agrees with all solar 
system tests of relativity provided F[I) appraches its linear form rapidly enough as / 



rises. It also reproduces much of the mass discrepancy phenomenology of spiral galaxies 
as encapsulated in Milgrom's MOND formula (Milgrom 1983a,b). Sanders (1986) has 
proposed a variant of AQUAL in which F{I) also turns linear in / for sufficiently low 
/. 

One problem with AQUAL and its variants is that the scalar field can propagate 
superluminally (Bekenstein 1987, 1988, 1992). This acausality is eliminated in the phase 
coupled gravity (PCG) theory (Bekenstein 1987, 1988, 1992). It has A = 1 and S = 
and instead of one scalar field, it has two, the usual ip and a second one A. The PCG 
dynamics comes from the action 

Si^,A = -\j g''^[AaAp+rj-^A^,a^,p+V{A^)]^d^x (8) 

where ?7 is a small coupling constant and y is a potential. For small rj the term in j[ 
involving i/j^a dominates the one involving A, a- Neglecting the latter term eliminates 
^ as a dynamical variable. Indeed, cxtrcmization of the action with respect to A 
establishes that ^4 is a function of g'-^^ip^a'^^jj^ so that to lowest order in ry, the PCG 
action reduces to AQUAL's, eq. (7). This suggests that PCG should be as succesful as 
AQUAL in explaining the phenomenology of the mass discrepancy. Several in-depth 
studies (Bekenstein 1987, Sanders 1988, 1989) have shown that this is the case, but have 
also uncovered potential problems with PCG in the solar system (Bekenstein 1987), and 
in galaxies (Sanders 1988a). 

From all the above examples we may write down the generic scalar action of interest: 

Si;,A = -J S{l,J,K:,A)V^d^x (9) 

Here £^ is a function, and X = g^f^ip,atp-,i3-, J = 9^^^ A,aA,p and /C = (7"'^^,q.^,^ are 
the three invariants that can be formed from first derivatives of A and (f). Although it 
would not change any of the discussion to come, ip is not included as an argument of 
£ because it would spoil the invariance of the theory under a shift of the zero oiip. At 
little cost we could generalize our action to include any number of ^-type fields. Note 
that we do not assume in what follows that the kinetic part of the scalar's action can 
be separated out, and is a quadratic form. We shall, however, assume that the function 
£^(X, J", /C, A) is such that the scalar field always bears positive energy density. In Sec. 5 
we translate this requirement into conditions on E and its derivatives. 

3. LIGHT BENDING IN GENERAL 

3.1. Equation for Light Rays 

In all modified gravity theories proposed so far as alternatives to dark matter, S = 
and the two metrics are taken as conformally related. Because Maxwell's equations are 
conformally invariant, this seems to say that the scalar field, because it enters only 
as a aconformal factor, has no infiuence on the propagation of light (Bekenstein 1992, 
Romatka 1992). Since the gravitational metric comes from Einstein-like equations, this 
would seem to say that all the ST theories with S = will give the same bending of light 



as GR. However, as mentioned in Sec. 1., observations of arcs suggest that the lensing 
in clusters of galaxies is much stronger than can be ascribed, via GR and formula (1), 
to the visible matter. This could mean that the modified gravity theories fail, and there 
is much dark matter in clusters. Alternatively, one can hope (Bekenstein 1992) that 
theories with B ^ could play a role in reconciling the large lensing with the idea of 
unconventional gravity by breaking the conformal relation between the metrics. Below 
we shall see that this hope is dashed because causality requires the "wrong sign" for B. 

The argument from conformality of the metrics to the absence of a light bending 
ascribable to the scalar field (Bekenstein 1992) leaves out an important factor. The 
scalar field, by virtue of its energy-momentum tensor, must make a contribution to 
gravitational metric g^p, and thus must lead to some extra light bending. The calcula- 
tions to follow are designed to take this extra bending into account. To accomplish this 
they are carried out at the fully relativistic level, rather than starting from results like 
formula (1) which are already stated in nonrelativistic terms. The results are surprising 
in that if the scalar field bears strictly positive energy, the extra light bending is nega- 
tive, i.e., if anything it decrements the general relativistic bending angle. As a result, 
no scalar-based gravity theory can make exactly the same predictions as GR for both 
dynamics of galaxies and clusters, and light bending. The opportunity to discriminate 
sharply between the two approaches thus arises. 

We shall discuss "light bending" of either electromagnetic or neutrino radiation. 
The Maxwell, massless Dirac and Weyl equations are all conformally invariant. Hence 
the physics of light bending may be discussed equally well by using the physical metric 
g(yp or the reduced metric 

9af3 = 9af3 + ^'^^a'^,p (10) 

where zu = L?'B{I) / A{I). It is important in all that follows that elementary considera- 
tions of signature and causality require that A{I) > while B{I) < (Bekenstein 1992, 
1993). Hence we may take w <0. 

For astrophysical systems the dimensions of the lensing system are large compared 
to the wavelengths of interest. Hence, we are only interested in geometric optics results. 
This means that our central question concerns the form of a null geodesic in the space- 
time defined by the reduced metric g^^p (we assume that "masses" like the neutrino 
mass or the plasma frequency in intergalactic space are small) . Therefore, the tracks of 
light rays through a gravitating cluster, x^^ = x^{^), where ^ is a parameter along the 
ray, must satisfy 

_ dx^dxl^ 

together with any restrictions arising from symetries and conservation laws. 

3.2. Bending in Spherical Symmetry 

It is sufficient for the points to be made to restrict attention to light bending by a 
static spherically symmetric configuration. As is well known, on symmetry grounds the 
metric of the corresponding spacetime, say g^^p, may always be put in the form 

ds^ = gapdx'^dx^ = -e^dt^ + e^dr'^ + {dO^ + sin O^dip^) (12) 



with V and A functions of r only. Similarly, we expect '4> = i(^{r) only. Hence eq. (11) 
takes the form 



+ +..^_J.O (13) 

where a prime denotes derivative with respect to r, and we have assumed that by angular 
momentum conservation, the path is confined to the equatorial plane. 

Conservation of energy (static configuration) tells us that (Misner, Thorne and 
Wheeler 1973) 

e^| = E (14) 

where E is the constant "energy" of the "photon" . We write "energy" because £"s value 
depends on the scaling of the parameter ^, which is arbitrary. Likewise, conservation of 
angular momentum (spherical symmetry) leads to 

.^|=. (15) 

with £ the "photon" 's "angular momentum" . £ is likewise changed by redefinition of ^. 
However, the ratio b = £/E is independent of such redefinition, and may be interpreted 
as the impact parameter of the ray from the center of the configuration. 

Substituting eqs. (14,15) into eq. (13) and replacing ^ by ^ as independent variable 
we cast the equation as 

Its first quadrature is 

J (e-£;2- £2/^2)1/2 ^2 ^'') 

This is the full description of the path of the ray in the form 9 — 9{r). The integral 
depends only on the ratio h = £/E, but not separately on E and £, thus realizing the 
independence on the parameter ^. 

How is the bending angle d determined ? Let ^ = stand for the direction of 
incidence of the ray being bent so that the lower limit of the integral in eq. (17) is 
r = oo. The integral from r = oo down to the turning point value, r = rturn where 
d9/dr = 0, with minus sign choice in the integral, is just half the rotation undergone by 
the radius vector as the photon comes in from infinity to closest approach. Were there 
no bending, the radius vector at closest approach would be at ^ = 7r/2. Hence, when 
there is bending, the net bending angle is given by 

Jn„ (e-nr/f))^ - 1)V2 r ^ 

with the TT representing the overall change of 9 when the path is perfectly straight. 

3.3. Bending in Linear Approximation 



Equation (18) is exact, but the integral cannot be evaluated analytically for the 
generic metric of interest. However, some approximations that can make the integration 
tractable are in order. Recall that A and v arise from a solution gfQ,^ of Einstein-like 
equations whose source is the Tq,^ of the matter in the system as well as of the scalar 
field generated by that matter. Because the system in question is weakly gravitating, we 
expect the metric g^^^ to be nearly fiat. Hence we may assume |A|, \v\ <^ 1. Likewise, 
the physical metric ^q,^ must be nearly flat by everyday experience in extragalactic 
astronomy. Of course, the flatness of the two metrics has to be evident in the same set 
of coordinates. Looking at eq. (3), we see that the conditions for common flatness of 
both metrics are that \'cu\ip''^ ^ 1 as well as A{I) exp{2if;) 1. From |A|, {wlip'"^ ^ 1 
we can obtain to linear order the following approximation to eq. (18): 



I..f(r2/62)n_,)_ 111/2 r 



It would, however, be premature to expand the integrand in v because an integral with 
(r^ — 52^3/2 ^Yie denominator diverges at the turning point. Appendix A shows how 
to handle this subtlety, and how to compute the integral to linear order in A, v and zu. 
The result is 



00 



r 



where an integration of the various quantities along a straight path with impact param- 
eter h is meant. Here z stands for linear distance. 

4. USING THE GRAVITATIONAL EQUATIONS 

4.1. Correction to Standard Theory Bending 

We now determine A, v from the Tq,^ distribution through the Einstein-like equa- 
tions. These are 

e-^(r-2 - r-^A') - r'^ = SttGT/ (21) 
e-\r-^v' + r-2) - r'^ = d,TTGT/ (22) 

Here 

Ta^ = ga^T^^, and T"/^, the energy-momentum tensor, is given by the variational 
derivative 

j.af3 = ^ '^('gm + S^) ^ ^^^^ 



The solution of eq. (21) is 



with 



-A 



l-2Gm{r)/r (24) 



rr 

m{r) = -Att j Tt^r'^dr (25) 

Jo 

Usually m(r) is interpreted as the gravitating mass up to radius r. Here, since Tt' 
is defined with respect to the gravitational metric and not to the physical one, m(r) 
is mass in unconventional local units which do not correspond to those measured by 



physical metersticks. Appendix B. develops the relation between the components of 
Tq,^ in the two frames of units. We find there that the factor that must be included in 
the integrand of eq. (25) to give the physical mass is very nearly unity for the situations 
of interest, so that m(r) is very nearly the physical mass. 

If we now expand e^^ in eq. (22) to 0(A), cancel a term, and neglect \v' jr in 
comparison with v' jr^ we get 

v> = \lr + 87rGrT/ (26) 
Substituting this and the linear approximation A = 2Gm{r)/r in (26) gives 

Note that Gm{r)/r'^ is just the magnitude of the formal Newtonian gravitational field 
g produced by all the mass energy interior to the radius r. The factor h/r is just the 
cosine of the angle between the radial direction and that normal to the ray at its nearest 
approach to the center. Hence the first integral in eq. (27) is identical to the standard 
prediction for -i?, eq. (1) (except for the fact that the scalar field's energy has to be 
included here in calculating g). The second integral gives the correction due to the 
scalar field which is not associated with its energy. 

We next recall that < by causality. Thus, even without knowing anything about 
the equation that determines V'j we can state categorically that 

/oo roo 
\g^\dz + AirGh j T/ dz (28) 
-oo J —oo 

We now proceed to break the integral over the stress T^^ into a contribution 6i} from 
the matter (stars, gas, dark) and another, Ai}, from the scalar field, and to show that 
the first is always negligible for the systems we have in mind. 

4.2. Evaluating the Stress Term 

The usual laws of energy-momentum conservation in special relativity must here 
take the covariant form 

T^,";,, = (29) 

where the semicolon signifies covariant differentiation with respect to the gravitational 
metric ^q,^. The law (29) follows from coordinate invariance of the action + -S'^ 
(Landau and Lifshitz 1975), or from the Bianch identities combined with the gravita- 
tional field equations derived by varying Sg + Sm + with respect to g^p. Now the r 
component of eq. (29) may be written in the form (Landau and Lifshitz 1975) 

(v^T/)' - \^g{dg^p/dr)T'^^ = (30) 

The assumed stationary and spherical symmetry forces ^q,^ and T'^f^ to be symmetric 
tensors. Thus eq. (30) takes the form 

(e^r^T/)' - ie^r2 \v'Tt' + A'T/ + 2(t/ + T^^)/r\ = (31) 



Spherical symmetry also demands Tq^ = T^^. Substituting this and noting that the 
terms involving A' cancel each other, we obtain the important result 

(e^/V^T/)' = ^e^/V^ pT/ + 4T//rj (32) 

As the radial component of the energy-momentum conservation equations, eq. (32) is 
the exact statement of hydrostatic equilibrium. 

Let us integrate this last equation over r in the interval [0,r]. Assuming that Tr'^ is 
bounded as befits a pressure, we see that the boundary term at r = vanishes. Hence 

T/(r) = ^ T^^I'^^v'Tt' + 4T//r)dr (33) 

In our weakly gravitating cluster we may replace the factors e^/^ by unity to get 

Tr'{r) ^^[^ r\u'Tt' + 4T//r)dr (34) 

which is our main working equation. 

4.3. Matter and Scalar Contributions 

There are two distinct contributions to Tr^ which enter ultimately into the expres- 
sion for light bending, eqs. (27,28). One belongs to the ordinary matter, and the other 
to the scalar field. First the ordinary matter (gas, stars, dark) contributes to and 
Tq. Because this matter has a velocity dispersion ~ f <^ 1, we expect its contribution 
to the stress \Tg \ to be of order v'^\Tt^\. To assess the magnitude of the first term in 
the integral in eq. (34), let us write down the physical line element that follows from 
eq. (12) in accordance with the transformation, eq. (3): 

It is evident from this that in the nonrelativistic limit (p = + ip plays the role of 
physical gravitational potential. The anomalous term ip here encapsulates the effects of 
unconventional gravity. Any theory of gravity whose task is to explain the mass discrep- 
ancy in terms of anomalous gravity must be such that ip' > (anomalous contributions 
strengthens the gravitational force). We may conclude from this that since v'^/r ~ 0', 
< u' < jr. Hence, according to eq. (34), the matter's contribution to Tr'ir) is 
kv^{p)r, where p denotes the matter's energy density, (— T/)ni, A; is a number of order 
unity and either sign, and the average used is defined by 



r"^ Jo 



rO dr (36) 

The scalar field makes equal contributions to and Tg in eqs. (33,34). For ac- 
cording to eqs. (9,23) the scalar's energy-momentum tensor is 

rj = -SgJ + 2{dS/dI)il;,aip,^ +2{dS/dJ)A,aA,^ +{dS/d]C){A,aip,^ +iP,aA,f^) 

(37) 



Because of the symmetries, neither ip nor A depend on t, 9 or (p. Therefore, 

V = 4 = t:^ = -8 (38) 

where it is understood that £ is to be evaluated using the solution to the field equations 
for if) and A. Our assumption (Sec. 2) of positive scalar field energy density obviously 
requires that —tu > 0, which in view of eq. (12) implies that r^* < and £ > 0, at least 
for a static situation. 

Because r^^ appears in eq. (34) multiplied by u' < v'^/r, the first term in the integral 
in eq. (34) is negligible compared to the second. Replacing in this last rj* by r^*, and 
recalling the definition of averaging, we see that the scalar 's contribution to 2r^(r) is 
{Tt*)r- We may summarize all the above in the equation 

T/(r) = kv^{p)r + {Tt% (39) 

where may vary from system to system, and even with r in accordance with the 
distribution of velocity and density, but should usually be of order unity. 

5. CONSEQUENCES OF POSITIVE ENERGY 

5.1. Scalar Field Decreases Bending 

To assess the size of the line integral over {p)r which enters into eq. (28) by virtue 
of eq. (39) we shall appeal to the Poisson equation for the Newtonian gravitational field 
gia generated by the matter. In cylindrical coordinates this is 

^ + ^ + ^ = -4-kGp (40) 
db b dz 

where we assume spherical symmetry to drop derivatives with respect to (f, and denote 
by b the coordinate in the direction of ±. Multiplying this equation by b and integrating 
over z gives 

+ / gidz = -AT,Gh pdz (41) 

where we have assumed that (7^ — ;> as |2;| — oo. Hence, taking into account that (p)^ 
is of the same order as p(&), we see that the contribution of the matter to the last term 
in eq. (28) is 

5d ~ -kv^ [b^ + 1^ j^^ gidz (42) 

Now, whatever the b dependence of the last integral, its logarithmic derivative in eq. (42) 
is unlikely to differ from itself by more than an order of magnitude. Since in a cluster 
is typically 10"^ km s"-*^ fa 3 x 10~"^c, it is clear that 5d may be neglected as compared 
with the first contribution to ■§ in eq. (28). (In a relativistic cluster, e.g. cluster of 
neutron stars at galactic center, v ^ 1 and would no longer be negligible). 

According to eqs. (28,36,39), the scalar field's contribution to the bending angle 
through the last term in eq. (28) is 

/oo J rr 
^ / rt*rdr (43) 
-00 ^ Jo 



It is already clear from the positivity of scalar field energy — r^* that A-i? < 0. As we 
shall see, this negative contribution dominates the positive one of the scalar field to the 
first term in {} as given by eq. (28) . We may transform the inner integral in eq. (43) by 
appealing to the Poisson equation for the formal Newtonian potential $s produced by 
the scalar's energy density. In spherical coordinates this is 

- (44) 

Multiplying through by r, integrating from r = to r, and noting that d^s/dr must be 
finite at r = 0, we find 

1^ r rMr = -"^^^ + (45) 

Of course, since comes from eq. (44), d^s{f)/dr = Gms(r)/r^, where ms(r) is the 
scalar's contribution to the total mass, c.f. eq. (25). Thus if we substitute eq. (45) in 
eq. (43) and that into eq. (28), and compare the last one with eq. (27), we see that the 
scalar's contribution to the first integral in eq. (28) is exactly cancelled out. We are 
thus left with 

^<2r|,^|d. + 26r5i(2^Mi,, (46) 

J—oo J—oo 

where gm is the formal Newtonian field of the matter alone. Now because the scalar's 
energy density — r^* > 0, it follows in the usual way from the Poisson equation (44) that 
grows with r. Thus the second integral in eq. (46) is strictly negative. Hence, 

/oo 
\gi\dz (47) 
-•oo 

This important result tells us that the scalar fields of the theory fail to augment the 
light bending ability of the matter, and may even weaken it. 

5.2. Divergent Gravitational Lenses ? 

What happens in the hypothetical case when the scalar field energy becomes the 
dominant part of the gravitational lense's mass ? We continue to assume that the 
configuration is nonrelativistic in the sense that A, etc. are small compared to unity. 
Then the same arguments as above lead to the cancellation of the scalar's contribution 
to the first integral in eq. (28). However, now we may neglect the matter's contributions 

and l^mM-^ as compared to the last term in eq. (45). The result is 

tf<ar?£MW,, (48) 

J-oo 

which makes it clear that the bending angle is negative. Such a scalar-energy dominated 
equilibrium configuration would thus behave as a divergent lense when bending light ! 

Under what conditions do equilibrium configurations dominated by scalar fields 
exist ? We may examine the question in the extreme case that matter is negligible 



without making the nonrelativistic approximation by returning to eq. (33), and replacing 
— > . By using eq. (38) we may recast eq. (33) into the form 

r/(r) = nr^e^/^yrM (49) 

Now the factor r^e'^/'^ should grow with r because the potential iy/2 is expected to 
increase with r (attractive gravitational force) , and because of the increasing factor r^. 
Even if over some range v were to decrease with r, it seems highly unlikely that e*^/^ 
would decrease faster than r^^. Thus positive energy density of the scalar field means 
here that at all radii Tr^ < 0. 

We now multiply eq. (49) by r^e^l'^ and differentiate with respect to r. The result 
may be put in the form 

(r/)' = {r\-l^)'e-l\-\n' - r/) (50) 
However, according to eq. (37) 



{d8/dI)^,l+{d8/dJ)A,l +{d8/diqA,r i^,r (51) 



We now argue that positive energy density of the scalar field with respect to any observer 
implies that the quantity in square brackets in eq. (51) is positive. If denotes the 
timelike 4-velocity of an observer, the scalar field will exhibit positive energy density 
to him if Tf^pW^U^ > 0. Since W^Up = —1, substituting eq. (37) into this condition we 
obtain the condition 



£ + 2 



{d£/dl){,p,a C/")2 + {d£/dJ){A,a f/")2 + id£/d}C)A,a t/"^',/? f^^l > (52) 



Now, we already know that £ > 0. Because t/" occurs in the quadratic form but not in 
£, positivity of this expression for any timelike t/" is guaranteed only if the quadratic 
form is positive definite. This requires 

d£/dl>0- dS/dJ>0 (53) 

and 

{dS/dJCf < 4{dS/dI){dS/dJ) (54) 

From these three conditions it follows that the quadratic form in the square brackets 
in eq. (51) is positive for any ip^r and A.,r- But then r^* — Tr'' < 0, so that it follows 
from eq. (50) that not only is Tr-'' negative for all r, but it also decreases outward. In 
particular, Tr'^{r) < — |rr^(0)| . This, of course, means that t^-'' can never reach zero, as 
would befit a finite configuration. This behavior also means, by eq. (51) that 

-n\r) > |t/(0) I + 2e-^ \{d£ / dl)tl;,l +{d£/dJ)A,l -^{dS/diqA^r > (55) 



which makes it clear that the energy density cannot vanish asymptotically as would 
be required for a configuration of finite mass. We conclude that a static spherically 
symmetric self-gravitating configuration of scalar field is untenable. It seems likely that 



the result would still hold for a configuration containing matter in which the scalar field 
dominates the mass of the matter. This makes the existence of divergent gravitational 
lenses very doubtful. 

6. CONCLUSIONS AND CAVEATS 

Comparing eqs. (1,47) we conclude that in a gravitational theory where gravity is 
mediated only by a metric and some scalar fields, the bending of light by a weakly gravi- 
tating system, like a galaxy or a cluster of galaxies, cannot exceed the bending predicted 
by GR for the mass of visible and hitherto undetected matter (but excluding the scalar's 
energy). This means that if gravity is properly described by a ST theory, and one uses 
the standard theory's formula (1) to interpret gravitational lensing observations, one 
can only underestimate the mass present in stars, gas and dark matter. 

The same conclusion obtains within GR if matter includes one or more scalar fields, 
i.e. Higgs fields. We obtain this case in our formalism by setting A{I) = 1, B{I) = and 
ip = 0. Then g^p — g^^j^. The A fields may stand for the scalar fields in question, and 
we may take over the results of our previous discussion: if gravity is described by GR, 
but matter includes various scalar fields, use of formula (1) to interpret gravitational 
lensing observations can only underestimate the mass of matter present not in scalar 
fields. 

Thus the observational claim that clusters of galaxies deflect light more strongly 
than would be expected from the observable matter contained by them, if it survives, 
cannot be interpreted in terms of scalar-tensor gravity without dark matter. Specifl- 
cally if follow-up observations eventually certify that the matter distribution inferred via 
standard theory from the lensing is very much like that determined from the dynamics 
of test objects or the temperature distribution of the X-ray emitting gas, then departures 
from standard gravity of ST form cannot play a role in the inner parts of clusters of 
galaxies which act as the lenses. 

A similar conclusion cannot be drawn for galactic gravitational lenses since lensing 
by single galaxies, if observed at all, is produced by the inner regions where, in all 
likelyhood, the visible matter is all there is. Neither is our last conclusion necessarily 
at variance with Milgrom's MOND scheme since the accelerations in the cores of rich 
clusters seem to be large compared to 10~^ c.g.s., the scale at which anomalous gravity 
effects set in according to MOND. In fact, what we have demonstrated is that ST 
unconventional gravity theories are irrelevant for understanding gravitational lensing 
by galaxy clusters. This does not exclude the possibilities that this type of theory may 
be of relevance for understanding the mass discrepancy in disk galaxies, or that the 
proper relativistic formulation of ideas like MOND involves more esoteric gravitational 
physics, e.g. modification of Newton's second law (Milgrom 1983a, 1993). 

In fact, it must be stressed that the very fact that a lensing mass, when deter- 
mined by standard theory, exceeds the directly detectable mass in stars or gas does not 
necessarily falsify unconventional gravity theories in general. There may be an as yet 
undetected component in clusters. A significant fraction, perhaps as large as 40%, of 
the mass in clusters, which was once considered to be dark, is now known to be in the 



form of hot gas. The matter in coohng flows must go somewhere, and such flows can 
deposit a substantial mass in the central regions of clusters. It would be interesting if 
clusters containing arcs all showed evidence for cooling flows. 

A system in which the classical dynamical mass, determined from virial methods or 
the distribution of hot gas, significantly exceeds the lensing mass as determined by GR, 
would be very problematic for the dark matter picture, but would be entirely consistent 
with unconventional ST gravity. This is not an entirely hypothetical situation. As 
mentioned in Sec. 1., the Tyson et. al measurements may be an example of this at the 
galactic scale. At cluster scales we may point out that what has been noticed in the sky 
are the striking examples of distant rich clusters containing luminous arcs or arclets; 
of equal importance may be the many examples of distant rich clusters without arcs or 
arclets (too low lensing masses). 

Our main results, e.g. eqs. (46,47) were obtained under the assumption that the 
scalar field configuration is truly static. There is thus a loophole in our conclusions. An 
equilibrium configuration containing a coherent complex scalar field oscillating harmon- 
ically in time, as in the boson star models (Liddle and Madsen 1992), is not covered by 
eq. (46). Neither is an equilibrium configuration in the PCG theory where scalar energy 
is nonnegligible, and where the cosmological expansion causes the field t/j to vary in time 
in approximately linear fashion (Sanders 1989). In future work we shall examine these 
situations. 

We both are grateful for the hospitality of the Aspen Center for Physics where this 
work was begun, and to Mordehai Milgrom for many conversations. J.D.B. thanks 
Rainer Romatka for comments and clarifications. 



APPENDIX A. 



Let us rewrite eq. (19) as 



d_ 
da 



roo 

/ (1 + A/2 + w^''^/2) [(r^b^) (1 -v)-a] ^l\-^dr 

^turn 



— TT 



a=l 



and then expand to first order in v. The results are 
with 



^1 



^2 



' 2 + A + wil)'^ dr 
b (r2/&2_ 1)1/2 7" 

Jb^a (r2/62-a)V2 



— TT 



(56) 



(57) 
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(59) 



We proceed to simphfy the expression for Two terms cancel out because of the 
integral 

Next, we transform the independent variable from r to 2; = ±(r2 — Ip')^/'^ which is just 
the Euclidean length along a ray whose impact parameter is 6, with zero length taken 
at the point of closest approach to the center, r = b. Thus 

X + wf^ 
-00 b^ + z"^ 



^1 



-dz 



(61) 



where we have extended the integral to 2; = —00 and divided by 2. 
We now turn to 'i?2- Integration by parts converts it to 



^2 = 



d_ 
b ] da 



roo 

lim iy(r2 -a62)V2 _ / j,//^2 _ ^^2a/2^, 



(62) 



a=l 



By expanding the square root under the limit in ab^/r^, we verify that after differ- 
entiation with respect to a, the corresponding term vanishes. Carrying out the a 
differentiation of the integral, setting a = 1, and passing from variable r to 2; we have 



^2 = 

Therefore, the full bending angle is 

b 



/oo 
-00 



1/ 



(62 + ^2)1/2 



dz 



(63) 



-'—00 \ 



A + wi/j 



/2 ^/ 

+ — 



dz 



(64) 



where again r = (b'^ + ^2)^/2. 



APPENDIX B. 

We begin by writing the tilde analog of eq. (23). 

In order to simplify the algebra we restrict ourselves to the case where A and B vary 
slowly, and neglect their derivatives. Since according to eq. (3) 

59ap/5~9ap = A-^e'^^ (66) 

we find by the chain rule that 

T"^ = (^/^)V2^-le-2V'Ta/3 (67) 

Of greater interest are the mixed components T^^ . Contracting both sides of eq. (67) 
with g^a and using the relation (3) and the definition w = B/Awe find that 

= (gf-g)^^^ (t/ + wT^^^,a^,^) (68) 

where indeces of T*^^ are lowered with g^p and those of T*^^ with g^p. 

The above result is general. For a static spherically symmetric situation the two 
metrics are given by eqs. (12,35) from which we infer that 

(^/^)l/2 = ^-2e-4^(l + e-A^^/2)-l/2 (gg) 

Further, because of the symmetries, tp,a = ipjr ^a^, and the mixed tensor T^^ is diagonal. 
Hence we find 

T//r/ = t//t/ = = A-2e-4V'(l + e-^ro^/^'2pl/2 ^^q^ 

and 

f//T/ = ^-2e-^^(l + e- W^)^/^ (71) 

Of course all the above results apply to the scalar field energy momentum tensor Ta^ 
alone. 

We may also use eq. (70) to reexpress eq. (25) for the mass function m(r) in terms 
of the physical energy density — T^*: 

m(r) = -47r A^e'^'f'{l + e-^wip'^)^/%*r^dr (72) 

Jo 

Since in the weak field situations under consideration Ae^V- ^ 1, A < 1 and < 1, 

it is clear that m(r) is very nearly the physical mass interior to radius r. 
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